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Abstract 

We consider the massive vector iV-component (Xip 4 )^ theory in Euclidian space and, using an ex- 
tended Matsubara formalism we perform a compactification on a ci-dimensional subspace, d < D. 
This allows us to treat jointly the effect of temperature and spatial confinement in the effective 
potential of the model, setting forth grounds for an analysis of phase transitions driven by temper- 
ature and spatial boundary. For d = 2, which corresponds to the heated system confined between 
two parallel planes (separation L), we obtain, in the large N limit at one-loop order, formulas for 
temperature- and boundary-dependent mass and coupling constant. The equation for the critical 
curve in the (3 x L plane is also derived. 



1 Introduction 



Questions concerning stability and existence of phase transitions are usually raised in the realm of 
field theories as an effect of temperature, but similar queries may also appear for confined systems, 
if we consider the behaviour of field theories as function of the spatial boundaries. One of the most 
striking examples of this kind of situation is the Abrikosov lattice of vortices occuring in type-II 
superconductors in high magnetic fields A recent account on aspects of finite-temperature 
field theory and on the behaviour of systems confined by spatial boundaries can be found in Refs. 
0, [3], [|, [5| . In another related domain of investigation, there are systems presenting defects created, 
for instance, in the process of crystal growth by some prepared circumstances (domain walls). At 
the level of effective field theories, in many situations this can be modeled by considering a Dirac 
fermionic field whose mass changes its sign as it crosses the defect. In other words, the domain 
wall plays the role of a critical boundary, separating two different states of the system (see for 
instance |7|, || and the references therein). In any case, the problem regarding to the behaviour 
of physical parameters in field theories (masses or coupling constants) on temperature as well as 
on spatial boundaries or inhomogeneities is still largely open and deserves further investigation. 
These aspects are addressed here by using a modified Matsubara scheme. 

For Euclidian field theories, imaginary time (temperature) and the spatial coordinates are on 
exactly the same footing, so that the Matsubara formalism can be also applied for the breaking 
of invariance along any of the spatial directions. Therefore, the existence of phase transitions 
induced either by spatial boundaries (separating different states of the system) or by temperature 
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or both can be analyzed in a general fashion. 

In order to advance with, in this note we present a general procedure to compactify a d- 
dimensional subspace for the massive vector iV-component (A<^ 4 )d theory in Euclidian space. 
This permits to extending to an arbitrary subspace results in the literature for finite temperature 
field theory 0, |3|, |J and for the behaviour of field theories in presence of spatial boundaries f5|. 
Our aim is to construct a general framework to study effects associated to simultaneous heating 
and spatial confinement, for a model that allows a non-perturbative approach. We consider the 
finite-temperature vector iV -component (A<£> 4 )i5 Euclidian theory at leading order in the system 
being submitted to the constraint of compactification of a d- dimensional subspace, d < D. After 
describing the general formalism, the renormalization is implemented for the interesting particular 
case with d = 2, i.e., the heated system at a temperature confined between two infinite parallel 
planes a distance L apart from one another. For this situation, we also derive the equation for the 
critical curve in the (3 x L plane. 

2 Compactification of a d- dimensional subspace 

We consider the model described by the Lagrangian density, 



in Euclidian D- dimensional space, where A is the coupling constant, m is the mass (see below) and 
summation over repeated "colour" indices a is assumed. To simplify the notation in the following 




(1) 
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we drop out the colour indices, summation over them being understood in field products. We 
consider the system in thermal equilibrium with a reservoir at temperature and confined 
to a (d — l)-dimensional spatial rectangular box of sides Lj, j = 2,3,...,d. We use Cartesian 
coordinates r = (xi, Xj, z), where z is a (D — d) -dimensional vector, with corresponding 
momentum k = (ki, kd, q), q being a (D — <i)-dimensional vector in momentum space. Then we 
are allowed to introduce a generalized Matsubara prescription, performing the following multiple 
replacements (compactification of a d -dimensional subspace), 

^ , g k ^ 

where it is understood that one of the axis (say, i — 1) corresponds to imaginary time (inverse 
temperature). A simpler situation is the heated system confined between two parallel planes 
a distance L apart from one another, normal to the x -axis. In this case we have Cartesian 
coordinates r = (r, x,z), where formally r correspond to the imaginary time axis and z is a 
(D — 2)-dimensional vector, with corresponding momenta k = (k T ,k x ,q), q being a (D — 2)- 
dimensional vector in momenta space. Then the multiple Matsubara replacements above reduce 
to, 

fdkr _ 1 +2? rdk x 1+2? 2n l7 r 2n 2 n 

Notice that compactification can be implemented in different ways, as for instance by impos- 
ing specific conditions on the fields at spatial boundaries. Here, we choose periodic boundary 
conditions in order to have a unified Matsubara prescription for both temperature and spatial 
dimensions. 



We start from the well known expression for the one-loop contribution to the zero-temperature 
effective potential 0, 



tfi(to) = E 



A0 



(4) 



(k 2 + m 2 ) s ' 

where m is the renormalized mass. For the Wick- ordered model, since in this case the tadpoles 
are suppressed, it is unnecessary to perform an explicit mass renormalization at the order at the 
one-loop approximation, the parameter m in Eq. flU) playing in this case the role of the physical 



mass. 

In the following, to deal with dimensionless quantities in the regularization procedures, we 
introduce parameters c 2 = m 2 /47r 2 /z 2 , (Lj/i) 2 = aj , g = (\/8tt 2 ), ((p /fi) = O , where y?o is the 
normalized vacuum expectation value of the field (the classical field) and \i is a mass scale. In terms 
of these parameters and performing the Matsubara replacements the one-loop contribution to 
the effective potential can be written in the form, 

oo / j^s +oo . ^D—dq 

Using a well-known dimensional regularization formula |§ to perform the integration over the 
(D — d) non-compactfied momentum variables, we obtain 

oo jj _ d 

Ux(<f>o, a u a*) = ^ D ^/aT 77T a~d^fi.D, d, s)g a 4^Af(s — ; a 1} a d ), (6) 

s=l 1 



where 



f(D, d, s) = „v-wLV^-r(s - (7) 



and 



A c d (u; a 1 ,...,a d ) 



+00 



(ain? H h a d r^ + c ; 

d 00 d 00 

i=l rii=l i<j=X rn,rij=X 

00 

+2 d ^ (ain? + • • ■ + a d r^ + c : 



1 

3i7 



2 1 2 1 2\-f 1 

+ a^n^ + c ) + 



„2\-z; 



(8) 



Next we can proceed generalizing to several dimensions the mode-sum regularization recipe de- 
scribed in Ref. JIU]. Using the identity, 



we get, 



where, 



A u V(y) Jo 



dt t u ~ l e-^\ 



A c d (u; ai,...,a d ) 



r(iz) Jo 



dt t v - x e 



l + 2]TT 1 (t,a,) + 



i=i 



+2 2 J2 T 2(t, <k, + + 2 d T d (t, at, a d ) 



(9) 



(10) 



ni=l 

Tj(t,ai,...,aj) = Tj-i(t,a\,...,aj-.{)Ti{t,aj) ,j = 2,...,d. 



(11) 
(12) 



Considering the property of functions Ti 



Ti(t,ai) = — + 4/ — 
u ; 2 V ad 



2 a,r 



(13) 
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where 



S(x) = J2 e~ n2x , (14) 

n=l 

we can notice that the surviving terms in Eq. fllOD are proportional to (a x ■ ■ -a^)^ 1 / 2 ). Therefore 
we find, 

rl 



A c d (v; ai, ...,a d ) 



d 
7T2 



Hdt £("-§)-l e -e 2 * 
JO 



d ,7T 2 



7T 2 ,7T 2 



d 7T 2 



• (15) 

i=l i<i=l fl * a j i=l 

Inserting in Eq. flloT) the explicit form of the function S(x) in Eq.flTil) and using the following 
representation for Bessel functions of the third kind, K v , 

2(a/b)%KJ2\/ab) = f° dx ^"V^"^ (16) 
Jo 

we obtain after some long but straightforward manipulations, 



A c d (u; ai,...,a d ) 



y/a,! ■ ^a~ d Y(p) 

d co 

+2EE 



2-l-i r ^--)(-) d -^ 

2 



in,=i V> 2 Lin 



)i- v K v _ § .(mL i n i ) + 



+2 d e 



rn 



ni,...,n d =l /! \J L^Tli + • • • + LjUq 



-^K v AmJLlnl + ... + L d n^ 



(17) 



Taking v = s — (D — d)/2in Eq . flT7|) we obtain from Eq.(|6] ) the one- loop correction to the effective 
potential in D dimensions with a compactified <i-dimensional subspace, 



Ui((f) ,ai, ...,a d 



s=l 



r-^Y(s--)C-) D - 2s 
2 



1 2 ' 



d oo 

i=ln,=l A" 
d oo 

+ 2 E E 



Kn_ s (mLirii) 



m 



i„ i ,n,=i> 2 A /L?nf + L 2 n 2 



-- S ^ £ _ (m,/L 2 n 2 + /. 2 //j) + ■ ■ 



+2 



d-1 



E 



i /iyi>? + • • • + i> 2 



18) 



with 



h{D,s) 



1 



:-i) 



s+1 



(19) 



2 D/2-s-l n D/2-2s s p( s ) • 

For d = 2, with Li = (3 and L 2 = we get the thermal and boundary dependent one-loop 
correction to the effective potential, 



T-f-*r( s --)(-) D - 2s 

2 n 



+2 E 



ni,n2- 



)f-: 



ifi2_ s (rW/3 2 n 2 + L 2 n|) 



(20) 



For = 1, Eq. flT8|) reduces to the cases of thermal or boundary dependences separately, studied in 
Refs. 0, Q and [BJ. In the next section we apply the formalism presented here to the simpler case 
d — 2: the heated system at a temperature confined between two parallel planes a distance L 
apart from one another. 
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3 Renormalization scheme 

We consider in the following the zero external-momenta four-point function, which is the basic 
object for our definition of the renormalized coupling constant. The four-point function at leading 
order in is given by the sum of all diagrams of the type depicted in Fig.l. This sum gives for 
the f3 and L -dependent four-point function at zero external momenta the formal expression, 

rS'(o,ftB = l 1 _ AS( A AftL) . (2D 

where H(D,/3,L) corresponds to the single bubble subdiagram in Fig.l. To obtain an expression 
for H(D,/3,L), we generalize one-loop results from finite temperature and confined field theory 
(H, M, |J) to our present situation. These results are obtained by the concurrent use of di- 
mensional and zeta -function analytic regularizations, to evaluate formally the integral over the 
continuous momenta and the summation over the Matsubara frequencies. We get sums of polar 
terms plus analytic corrections depending on the compactified variables. Renormalized quantities 
are obtained by subtraction of the divergent (polar) terms appearing in the quantities obtained by 
application of the modified Feynman rules (Matsubara prescriptions) and dimensional regulariza- 
tion formulas. These polar terms are proportional to T-functions having the dimension D in the 
argument and correspond to the introduction of counter terms in the original Lagrangian density. 
In order to have a coherent procedure in any dimension, these subtractions should be performed 
even in the case of odd dimension D, where no poles of T- functions are present. 

Note that since we are using dimensional regularization techniques, there is implicitly in the 
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above formulas a factor /i 4 ~ D in the definition of the coupling constant. In what follows we 
make explicit this factor, the symbol A standing for the dimensionless coupling parameter (which 
coincides with the physical coupling constant in D = 4). To proceed we use the renormalization 
conditions, 



and 



ff 2 



-U 1 (D,f3,L)\ (l , 0=0 = m 2 n 2 



(22) 



(23) 



from which we deduce that formally the single bubble function H(D,f3,L) is obtained from the 



coefficient of the fourth power of the field (s = 2) in Eq. (20). Of course such coefficient is in 
general ultraviolet divergent and a renormalization procedure is needed. Then using Eqs. (|23"D 
and (|20"D we can write Y>(D,/3,L) in the form, 



Z(D,(3,L)=H(D)-G(D,P,L), 



(24) 



where the f3 and L dependent contribution G(D, (3, L) comes from the second term between brack- 
ets in Eq. (|20|) , that is, 

3 /i 4 - D 



G(D,f3,L) 



-2 E 



m 



ni,n 2 =l \l p 2 n\ + L 2 n\ 



(25) 



and H(D), a polar parcel coming from the first term between brackets in Eq.(^), is 

H(D)KT(2-^)(m/v) D -\ 



(26) 



We see from Eq.([26|) that for even dimensions D > 4, H{D) is divergent, due to the pole of 
the T-function. Accordingly this term must be subtracted to give the renormalized single bubble 
function T,r(D, j3, L). We get simply, 

Z R (D,(3,L) = -G(D,(3,L). (27) 

As mentioned before, for sake of uniformity the term H(D) is also subtracted in the case of odd 
dimension D, where no poles of T -functions are present. In these cases we perform a finite 
renormalization. ^From the properties of Bessel functions it can be seen from Eq. (|25|) that for any 
dimension D, G(D,j3,L) satisfies the conditions, 

lim G(D,P,L) = 0, lim G(D,P,L) -> oo, (28) 

and G(D, (3, L) > for any values of D, (3 and L. 

Let us define the (3 and L dependent renormalized coupling constant \r(D, (3, L) at the leading 
order in 1/N as, 

r£'(o. p,l) . 1a k( aM = jf r^i^i) (») 

and the renormalized zero-temperature coupling constant in the absence of boundaries as, 

^Xr(D)= lim rg>,/?,L). (30) 

iV P,L— >oo 



From Eqs.flSTp, (P^), (PH) and (|2T|) we conclude that Xr(D) = X. In other words we have 



done a choice of renormalization scheme such that the constant A introduced in the Lagrangian 
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corresponds to the physical coupling constant in free space, at zero temperature. From Eqs.(p9|) 
and (|27D we obtain the (3 and L dependent renormalized coupling constant, 

A«(A/^ 1 + AG( A W . (31) 

The above procedure corresponds, on perturbative grounds, to sum up all the chains of bubble 
graphs in Fig.l. It is just the resummation of all the perturbative contributions including the 
counter terms from the chains of bubbles and the subtraction of the divergent (polar) parts, written 
in compact form. Again, these subtractions should be performed even for odd dimensions. 

We have used above a modified minimal subtraction scheme where the mass and coupling 
constant counter terms are poles at the physical values of s. The multiple and double sums in 
Eqs. (|8|) to ( |2l]| ) define generalized inhomogeneous Epstein-Hurwitz zeta-functions, in such a way 
that the (3 and L dependent correction to the coupling constant is proportional to the regular 
part of the analytical extension of these inhomogeneous Epstein-Hurwitz zeta-functions in the 
neighbourhood of the pole at s = 2. For the non- Wick- ordered model the same argument applies 
to the renormalized mass for s — 1, and the f3 and L dependent renormalized mass at one-loop 
approximation is given by, 



oo oo 
Ln=l U P - n=l UL 

m 



+2 £ ( 7= 2 r2 ^Kn^my/Pnl + Vnl 
ru,n a =i y P ' n{ + L l n{ 



(32) 



In dimension D = 3, from the relation 11 



= V^ 6 "*' (33) 
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after summing up geometric series, we obtain the following semi closed expressions 

3/i 



G(P,L) 



8nm 



[p m P _ I)" 1 _|_ [p mL — l)- 1 _|_ 2 y~] e -™^/ I3 2 n\+L 2 n. 



ni,ri2=l 



(34) 



log(l - e- m/3 ) log(l - e~ 



(3 



+ 2 



oo 

E 



i^Jp 2 n\+L 2 n 



4 Mass behaviour and critical curve 



(35) 



As far as the mass behaviour is concerned, if we start in the ordered phase with a negative squared 
mass, the model exhibits spontaneous symmetry breaking of the O(N) symmetry to 0(N — 1), 
but for sufficiently small critical values of j3 and L the symmetry is restored. We can define the 
critical curve C((3 C ,L C ) as the curve in the j3 x L plane for which the inverse squared correlation 
length , £~ 2 (/3, L, (p ), vanishes in the gap equation, 



ml + Xifl 

MN + 2) 
2N(3L 



E 



ni,rt2=— oo 



d 



D-2, 



(2ir) D - 2 q 2 + (2Ha) 2 



;^) 2 + e- 2 (/?,4^ ) 



(36) 



where ipo is the normalized vacuum expectation value of the field (different from zero in the ordered 
phase). 

Close to the critical curve, <po vanishes and the thermal gap equation at one-loop order reduces, 
in the large N limit, to Eq. (|32|) . If we limit ourselves to the neighbourhood of criticality, m 2 w 0, 
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we may use an asymptotic formula for small values of the argument of Bessel functions, 



KJz) 



;*~o) 



(37) 



and Eq.([32f) becomes, 



m 2 (/3, L) « ml + " v ',, , H — - I ) 



2 



7T 



,2J/2 



09 



2--D , T 2-D 



L 2 - U )C(D - 2) + 2£ 2 ((L> - 2)/2; ^, L ) 



(38) 



where £(.D — 2) is the Riemann zeta-function, 



c(d-2) = y: 



n=l 



n 



D-2 



(39) 



and E 2 ((D — 2)/2; f3 2 , L 2 ) is the generalized Epstein-Hurwitz zeta -function |L2 

E 2 ((D -2)/2; f3 2 ,L 2 ) 



E 

ni,n2=l 



1 



(/3 2 nf + L 2 n|)^ 



(40) 



The Riemann zeta-function ((D — 2) has an analytical extention to the whole complex Z)-plane, 
having an unique simple pole (of residue 1) at D = 3. On the other hand, from the properties of 
generalized Epstein zeta -functions, in particular recurrence relations, we can write the generalized 
zeta-f unction fl4"0l) in terms of the Riemann zeta-f unction. As it is done in Ref . [|12| , whatever the 
sum one chooses to perform firstly, the manifest (3 <-> L symmetry of Eq. 



is lost. In order to 

preserve this symmetry, we adopt here a symmetrized summation which leads to 
E 2 ((D-2)/2;f3 2 ,L 2 ) = ~ (/3 2 ~ D + L 2 ~ D ) ((D - 2) 



4r( 



D-2 



+ 



CP -3) 



D-2 



-W 2 {D-2;P,L) 



(41) 
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where 

W 2 (D-2;3,L) = -— ]T (7r-n 1 n 2 )^ J ftT^3(27r-n 1 n 2 ) 

P ni,n 2 =l P P 

1 00 3 d 3 
+ Td^ £ (^j-n 1 n 2 )^KD-3(27r^nin 2 ) , (42) 

rai,n.2=l 

from which we conclude that E 2 ((D — 2)/2; 3 2 , L 2 ) has only two simple poles at D = 4 and 
D = 3. The solution of Eq.(^) for m(/9, L) = and 7Bq < defines the critical curve C(/3 C , L c ) in 
Euclidean space dimension D (D 7^ 4 and Z) 7^ 3), 



4Xfi A - D ^.D 



m o + T^TT^y ~ !) [(#" + L c KP - 2) + 2£ 2 ((D - 2)/2; L e )j = 0. (43) 

For D = 4 the generalized ,zeta-function has a pole and for D = 3 both the Riemann zeta- 
function and the function Z 2 have poles. We can not obtain a critical curve in dimensions D = 4 
and D = 3 by a limiting procedure from Eq.([|3|). For D = 4, which corresponds to the physical 
interesting situation of the heated system confined between two palallel planes embedded in a 
3-dimensional Euclidean space, Eqs.(^) and (f43"|) become meaningless. To obtain a critical curve 
in D = 4, a regularization procedure can be done as follows: we perform an analytic continuation 
of the zeta-function to values of the argument z < 1, by means of the formula, 

Eq . (|44]) defines a meromorphic function having only one simple pole at z — 1. Next, remembering 
the formula, 



lim 

Z->1 



CO*)" ' 



1-2 
14 



C , (45) 



where C is the Euler-Mascheroni constant, and using Eq. ([41]) , we define the renormalized mass 



lim 7ftn + 
D^4_ u 



1 4A 




(46) 



D-4n 2 (3L 



in terms of which we obtain the critical curve in dimension D = 4, 




7rC 



+ A v ^W 2 (2;f3 C} L c )=0 1 



(47) 



where we have used C(2) = 7r 2 /6 and T(l/2) = y^. 

5 Concluding remarks 

In this paper we have presented a generalization of the Matsubara formalism to account for 
simultaneous compactification of imaginary-time and spatial coordinates. The method is applied 
to the massive vector N-component (A<^ 4 )d theory, in Euclidian space, compactified in a d - 
dimensional subspace (d < D), including temperature. We have calculated the effective potential 
of the model from which, in the particular case d = 2 , equations for the (3- and L-dependent 
mass and coupling constant have been derived. In the absence of Wick-ordering, the renormalized 
mass can not be taken as the coefficient of the term ip a fa in the Lagrangian, and we must take 
a j3- and L- corrected renormalized mass, m(f3,L). To obtain the (3- and L- dependent coupling 
constant, the mass parameter m should be replaced by m(/3, L) in Eqs( [H]) and (^) and the 
resulting system of equations should be solved. Under these conditions exact closed expressions 
are almost impossible to obtain, since it would require a procedure equivalent to solving exactly 
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the Dyson-Schwinger equations. For D = 4, an equation for the critical curve in the (3 x L plane 
is obtained. An application of the formalism presented here to perform a detailed analytical and 
numerical study of the behaviour of the coupling constant and the mass on temperature and on 
spatial boundaries is in progress and will be presented elsewhere. 
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Figure captions 

Fig.l - Typical diagram contributing to the four-point function at leading order in j^. To each 
vertex there is a factor and for each single bubble a colour circulation factor N. 
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